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Abstract 

We extend to the first 300 Riemann zeros, tlie form of analysis reported by us in Canad. J. 



Phys. 85, 1 (2007) (also arXiv:math-ph/0606005 ), in which the largest study had involved the 



first 75 zeros. Again, we model the nonsmooth fluctuating part of the Wu-Sprung potential, 



i-^ . which reproduces the Riemann zeros, by the alternating-sign sine series fractal of Berry and 

Ph' 

r~| ! Lewis A{x,^). Setting the fractal dimension equal to |, we estimate the frequency parameter (7), 

» . 

^ ' plus an overall scaling parameter [a] introduced. We search for that pair of parameters (7, a) 

which minimizes the least-squares fit of the lowest 300 eigenvalues — obtained by solving the 
one-dimensional stationary (non-fractal) Schrodinger equation with the trial potential (smooth 
' plus nonsmooth parts) — to the first 300 Riemann zeros. We randomly sample values within the 

o ' 

ff^ ■ rectangle < o" < 3, < 7 < 25. The fits obtained are compared to those gotten by using simply 



the smooth part of the Wu-Sprung potential without any fractal supplementation. Some limited 



o 

I improvement is again found. There are two (primary and secondary) quite distinct subdomains. 



in which the values giving improvements in fit are concentrated. 
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I. INTRODUCTION 



Our objective here is the discernment of structural patterns within the nontrivial Rie- 
mann zeros (that is, the real parts > of the zeros of the Riemann zeta function C{s), 
which conjecturally all take the form ^ + aki)- More specifically, following our work in 







l| (also arXiv:math-ph/0606005), we seek to model the fractal component of the potential 



which reproduces the a^s. The (dominant) smooth part of such a potential can be written 
implicitly in the form, 

.WV')^i(Vl^..^.VF>„f±^). (1, 



(Here Vq = S.IOOTStt ^ 9.74123.) This formula (HD was derived by Wu and Sprung [2j by 
solving a form of Abel's integral equation obtained by the application of the semiclassical 
approximation to the formula, 

n 

for the smooth part of the number of Riemann zeros less than E (cf. 131). 

n 

In pj, we modeled the fractal component by the alternating-sign sine series fractal (a 
particular case of a deterministic Weierstrass-Mandelbrot fractal function) of Berry and 
Lewis 1^, eq. (5)], 

A(x,7) = S-_^ ^~y^^J"^^ (l<d<2,l<7). (3) 

Here, d is the fractal dimension, which — following the box-counting argument of Wu and 



Sprung (cf. [5[]) — we took to be |. (van Zyl and Hutchinson also came to the same 
conclusion regarding the fractal dimension, independently of which of two distinct inversion 
techniques they employed We have, in this d = ^ Berry-Lewis context, a specific case, 

A(7x,7) = -7^A(x,7), (4) 

of the "affine scaling law" {4, eq. (3)]. _ 

In the most extensive of the three (n = 25,50,75) analyses reported in [l], we found 
that values of the frequency parameter 7 in the vicinity of 3 (and, secondarily, 9) tended to 
provide the best fits to the Riemann zeros (a^'s) themselves (ll. Fig. 16]. (The sum of the 
inverted Wu-Sprung potential ([1]) plus the alternating-sign sine series ([3]), with 7 fixed at 



various random values, was used in these analyses as an adjusted potential, which served 
as input to the time-independent one-dimensional Schrodinger equation. The eigenvalues 
obtained are, then, regarded as fits to the corresponding Riemann zeros.) 

Here, we continue along these same investigative lines, quadrupoling the sample size of 
zeros from 75 to 300. We achieve this substantial increase mainly through reducing — from 
5,625 to 600 — the basis size employed in the Arnoldi method for finding eigenvalues of 
sparse matrices as well as sampling 6,500 points, not 10,000, in the Mathematica program 
supplied to us by M. Trott. (We also truncate the series ([3]) at m = 40, while in [l|] we used 
a cutoff of m = 30 — which, of course, does not reduce computational time.) Despite these 
reductions, our program, when we employ just the smooth (unsupplemented) Wu-Sprung 
potential ([T]), yielded what we deemed to be a suitably acceptable prediction of 542.023 for 
the (highest) 300-th Riemann zero, that is, 541.847. ( It would certainly be of interest to 
evaluate the sensitivity of our results to the various computer program parameters employed 
and adjust them accordingly, but it would seem to be overwhelmingly time-consuming to 
do so in substantial detail. We also tried to adopt our Mathematica program to analyze 
the first 500 Riemann zeros, but we could not adjust the parameters to provide us with an 
acceptable fit, as measured by comparing the 500-th eigenvalue with the 500-th Riemann 
zero.) 

While, in the n = 75 analysis in [l|, we randomly sampled values of 7 G [0, lOj^here we 
expand the interval to 7 G [0,25], and as in the first of the analyses (n = 25) in p|, allow 
the series ([3]) to be scaled by an overall factor of a, we take now to lie in the interval [0,3]. 
(We arrived at this rectangular range of parameters from which to extensively sample, after 
a number of preliminary analyses.) 

In Fig. [1] we show the results obtained by randomly sampling points in the indicated 
base two-dimensional rectangle and determining how well the so-supplemented smooth Wu- 
Sprung potential fits the first 300 Riemann zeros. Only those 213 points (of the 3,330 
sampled) for which the fit is improved for all three measures Skier, •y) = S^^^^IcJa; — afcp, j = 
1,|,2 are displayed. (The a^'s denote, say, the eigenvalues, and the a^'s the Riemann 
zeros.) These measures were, respectively, 116.265, 86.3596 and 68.1307 when no fractal 
component was added, and only the smooth potential employed. So, Fig. [1] contains only 
points for which all three indicators are less than these three values. All the improvements 
are rather minor in character — as in [l| — with none exceeding one-percent in terms of the 
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FIG. 1: Those 100 pairs (cr, 7) which give improved alternating-sign sine series ([3]) fits (for all 
k = 1, |,2) to the first 300 Riemann zeros. 20,000 pairs were randomly, uniformly sampled over 
the base rectangle (0 <cr< 3, 0<7< 25). The third axis 5*2 (cr, 7) gives the associated least-square 
fit (j = 2) — smaller values, of course, denoting better fits to the zeros. The minimum value of 
67.562 is attained at (10.4444, 1.50996). 

criterion values ^^(cr, 7). We observe that for most of these points yielding improvements 
in fit, we have 7 > 10, a region that was outside the range of consideration in The 
minimum values using the j = I and J = | criteria, 115.451 and 85.6691, were attained at 
the same point 7 = 10.3954, a = 1.55535, while using the j = 2 (least-squares) criterion, the 
minimizing point, with the value 67.562, was located at 7 = 10.4444,0" = 1.50996. (For all 
of the 213 improving values, we had a < 1.85799.) 

In Figs. [2] and [HI we magnify the regions in Fig. [H in which the improvements in fit are 
concentrated. 

We also investigated the possibility of modifying our analyses by adding a prefactor of 
7~t to ([3]) to enhance periodicity, following certain suggestions of Tricot ^, sec. 12.11], but 
this led to convergence problems, which we did not know immediately how to overcome. Of 
course, it would be of interest to employ still other basic (few-parameter) fractal models in 
addition to ([3]) in similar analyses of the fractal component of the Wu-Sprung potential for 
the Riemann zeros, possibly thereby achieving greater improvements in fits to the zeros than 
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FIG. 2: Magnification of primary domain of fit-improving points in Fig. [T] 




FIG. 3: Magnification of secondary domain of fit-improving points in Fig. [T] 

we have so far been able to obtain. (In this regard, the companion cosine series of Berry 
and Lewis {4, eq. (4)] is of potential interest, although unhke ([3]), it is never negative, so it 
would seem necessary to add some negative quantity to it.) 

It would also be of interest to directly model the fractal component remaining when the 
smooth Wu-Sprung potential is fitted to the Riemann zeros (cf. 0, inset. Fig. 2] j^. Fig. 3 
(a)-(c)]), rather than having — as we have done here and in [l| — inputting each candidate 
for the fractal component into the time-independent one-dimensional Schrodinger equation. 
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